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Foreword 


This  report  describes  work  that  was  carried  out  in  the  Thermosciences  Division 

at  NASA  Ames  Research  for  the  period  1990-1992.  The  work  is  the  combined 
effort  of  the  Computational  Chemistry  and  Aerothermodynamics  Branches. 

The  principal  goal  of  the  work  was  the  prediction  of  radiation  signatures  from 

shock  layers  from  ascending  vehicles.  Some  work  was  also  directed  at 
understanding  the  chemistry  in  the  plume  and  plume  induced  separation.  The 

research  in  the  two  branches  was  complementary  in  that  both  accurate 
molecular  properties  and  the  flow  characteristics  must  be  included  in  the 
codes  that  attempt  to  produce  the  signatures  observed  on  actual  flight 
experiments.  Nevertheless,  the  work  in  the  two  branches  is  sufficiently 

distinct  that  the  contributions  from  the  two  are  described  separately  in  the 
final  report. 
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Computational  Chemistry  Branch 
A.  Statement  of  the  Problem 

The  ultimate  goal  of  the  research  is  to  be  able  to  predict  the  radiation 
signature  of  ascending  vehicles.  This  is  essential  in  order  to  intercept  missiles 
during  their  boost-phase.  Our  recent  work  is  directed  at  demonstrating  that 
detection  of  the  hard  body  in  the  ultraviolet  region  where  the  Earth’s 
background  radiation  is  small  is  a  viable  concept.  There  are  two  actual  flight 
experiments  that  have  provided  very  accurate  spectral  data  in  the  3-5  km/s 
velocity  range.  In  collaboration  with  the  aerothermodynamics  branch  at 
Ames  as  well  as  other  researchers  in  this  program  we  have  worked  on 
understanding  the  physics  of  shock  heated  air  and  in  the  development  of  codes 
to  predict  the  radiation  signature  as  a  function  of  altitude  and  speed.  This 
section  of  the  final  report  describes  the  work  that  was  done  in  the 
Computational  Chemistry  Branch. 

The  molecular  data  that  is  needed  to  predict  the  bow-shock  radiation  and  to 
understand  the  chemistry  and  emission  from  the  plume  can  be  divided  into 
five  major  classes:  (1)  radiative  transition  probabilities  for  molecular 
transitions  that  contribute  to  the  signatures;  (2)  transport  properties  such  as 
thermal  diffusion,  thermal  conductivity  and  viscosity  needed  to  describe  the 
flow  in  the  shock;  (3)  thermodynamic  properties  such  as  heats  of  formation  to 
understand  the  chemistry  in  the  plume;  (4)  cross  sections  for  electron 
molecule  interactions;  and  (5)  reaction  rates  for  both  equilibrium  and  non¬ 
equilibrium  reactions. 
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One  of  the  main  objectives  of  our  research  program  is  to  compute  radiative 
intensity  factors  for  air  species  that  contribute  to  the  radiation  signatures  of 
hard  bodies  and  exhaust  plumes.  The  NO  molecule  is  one  of  the  key  species,  as 
it  is  formed  in  significant  concentrations  in  shock  layers,  for  example,  from 
the  reaction  of  oxygen  atoms  with  N2-  During  the  three  year  period  of  this 
proposal  we  have  determined  accurate  transition  probabilities  for  the  beta, 

ogawa,  and  epsilon  systems  of  NO.  The  beta  system  (B^ri'-X^n)  was  particularly 
challenging  theoretically,  because  the  upper  state  changes  from  Rydberg  to 
valence  character  in  the  Franck-Condon  region  of  the  transition. 
Nevertheless,  the  radiative  lifetimes  as  a  function  of  vibrational  level  are  in 
good  agreement  with  recent  laser  induced  fluorescence  measurements.  The 
theoretical  results  for  the  NO  ogawa  system  (b^Z‘--a^n)  do  not  agree  as  well 
with  experiment.  However,  we  believe  the  theoretical  values  are  the  most 

accurate  available  for  this  band  system. 

More  recently  we  have  been  computing  accurate  radiative  intensity  factors 
for  the  epsilon  system  of  NO.  The  epsilon  bands  are  a  prominent  constituent  of 
shock-heated  air  in  the  spectral  region  between  180  and  210  nanometers.  Our 
theoretical  radiative  lifetime  of  19  nanoseconds  for  the  low-lying  vibrational 
levels  is  in  good  agreement  with  experiment.  Thus  our  Einstein  coefficients 
for  this  band  system  should  represent  a  substantial  improvement  over  existing 

values.  These  data  are  being  incorporated  into  the  NEQA1R  code  for  computing 

the  radiation  from  shock-heated  air. 

Emission  from  the  Cameron  band  system  of  CO  is  observed  in  the  plumes  of 

rockets.  To  gain  insight  into  the  possible  mechanism  for  this  long-lived 

chemiluminescence,  we  have  studied  the  quintet  states  of  CO  that  are  populated 
from  recombination  of  ground  state  atoms.  Of  all  the  quintet  states  dissociating 
to  ground  state  atoms,  only  two  are  found  to  be  significantly  bound.  This  work 
resulted  in  the  assignment  of  two  new  band  systems  of  CO,  and  it  was  shown 

that  the  emission  observed  between  about  360  and  520  nanometers  is  due  to  the 

a  ^n--a^E+  and  a  ^n--d^A  transitions.  Thus  it  is  possible  that  the  quintet 
states,  populated  through  recombination  of  ground  state  atoms,  serve  as 
precursor  states  for  the  formation  of  the  upper  state  of  the  Cameron  system. 
However,  the  lifetime  of  this  process  is  too  short  to  explain  the  long-lived 
chemiluminescence  from  the  Cameron  system. 

We  computed  the  thermodynamic  properties  of  various  A1-C1-0  radical  species 
that  are  expected  to  be  present  in  rocket  plumes.  Specifically  we  determined 
the  thermochemical  heats  of  formation  of  plume  species  comprised  of  one  or 
two  aluminum  atoms  and  various  numbers  of  chlorine  and  oxygen  atoms. 

Work  is  continuing  on  refining  the  geometry  of  the  most  stable  isomer  of 
AI2O2  with  an  extensive  treatment  of  electron  correlation. 
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Our  laboratory  has  been  a  leader  in  the  computation  of  accurate  transport 
properties.  During  the  last  three  years  we  considered  the  interaction  of 
nitrogen  and  oxygen  atoms  and  ions.  The  calculations  are  based  on  ab  initio 
potentials  and  a  semiclassical  approximation  to  the  collision  integrals.  This 
latter  approximation  has  been  shown  to  be  valid  at  room  temperature  and 
above  based  on  comparisons  with  an  accurate  quantum  treatment  of  the 
scattering  problem.  These  atom-atom  and  atom-ion  transport  properties  such 
as  thermal  diffusion  and  viscosity  are  thought  to  be  accurate  to  5%. 
Calculations  are  in  progress  to  determine  accurate  atom-molecule  and 
molecule-molecule  collision  integrals  for  important  air  species. 

Electron-molecule  collisions  become  an  important  mechanism  for  exciting 
molecules  at  high  velocities.  These  collisions  are  particularly  efficient  at 
scrambling  the  populations  in  the  excited  states.  Thus  to  predict  the  emission 
from  the  shock  layer,  it  is  necessary  to  determine  accurate  cross  sections  for 
electron-molecule  collisions  between  the  excited  states  of  important  air 
species.  In  this  proposal  period  we  determined  cross  sections  for  some  excited 
states  of  N2- 
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B  Summary  of  Most  Important  Results 

1.  Accurate  transition  probabilites  and  branching  ratios  for  the  beta,  ogawa, 
and  epsilon  band  systems  of  NO  were  determined.  These  are  the  most  accurate 
presently  available  and  these  data  have  been  incorporated  into  the  NEQA1R 
codes  used  to  predict  the  radiation  from  shocks. 

2.  Assigned  several  new  band  systems  of  the  CO  molecule  involving  transitions 
between  the  quintet  and  triplet  states.  Demonstrated  that  it  is  unlikely  that  the 
quintet  states  of  CO  could  provide  an  excitation  mechanism  for  the  upper  state 
of  the  Cameron  band  system  observed  prominantly  in  plumes. 

3.  Accurate  transport  properties  were  computed  for  atom-atom  and  atom-ion 
interactions  involving  nitrogen  and  oxygen  atoms.  Collision  integrals  and 
transport  properties  were  reported  in  a  series  of  papers  for  a  wide  range  of 
temperatures. 

4.  Electron  collision  cross  sections  were  computed  for  some  excited  states  of  N2. 
which  are  important  in  determining  the  state  populations  in  the  shock  layer, 
especially  at  the  higher  velocities.  We  also  computed  vibrational  excitation 
probabilites  for  N2-N2  collisions,  which  is  a  key  mechanism  for  accomodating 
the  energy  in  the  shock  layer. 
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Aerothermodynamics  Branch 

A.  Statement  of  the  Problem 

The  ultimate  goal  of  the  research  within  the  Aerothermodynamics  Branch  is 
the  same  as  that  for  Computational  Chemistry  Branch.  The  specific  goal  of 
Aerothermodynamics  Branch  is  to  develop  a  computer  code  that  numerically 
reproduces  the  flight  data.  This  consists  of:  (1)  interpreting  the  flight  data.  (2) 
developing  the  theoretical  model  that  describes  the  thermochemical 
phenomena  occurring  in  the  flow,  (3)  developing  the  CFD  code  that 
accommodates  the  thermochemical  model,  (4)  run  the  code  for  the  flight 
conditions,  and  compare  the  calculated  results  with  the  flight  data. 

Aerothermodynamics  Branch  has  grants  to  both  Stanford  University,  with 
Prof.  Dean  Chapman  as  the  principal  investigator,  and  North  Carolina  State 
University,  with  Prof.  Graham  Candler  as  the  principal  investigator,  to  carry 
out  portions  of  the  task.  These  two  universities  have  had  other  professors,  such 
as  Prof.  Robert  MacCormack,  and  several  graduate  students  to  assist  the 
principal  investigators.  Ames  staff  collaborated  with  these  university 
personnel.  In  addition,  Ames  staff  collaborated  also  with  Institute  for  Defense 
Analysis  (IDA),  with  Dr.  Debrah  Levin  as  the  principal  investigator,  which  has 
been  funded  independently  of  Ames. 

The  Bow  Shock  Ultraviolet  Flight  1  was  flown  in  April  1990,  at  a  flight  of  about 
3.5  km/s.  The  spectrometers  and  monochromators  located  in  the  nose  region  of 
the  vehicle  measured  a  substantial  radiation  in  the  ultraviolet  region,  which 
consists  mostly  of  nitric  oxide  (NO)  band  systems.  The  Flight  2  was  flown  in 
February  1991  at  a  flight  speed  of  about  5.1  km/s.  In  the  second  flight,  the 
radiation  from  the  plume  region  was  determined  as  well  as  the  nose  region. 

One  radiometer  viewed  the  plume  approximately  along  the  centerline,  while 
the  other  viewed  far-field,  that  is,  at  25  degrees  off  the  centerline.  The 
intensity  of  NO  radiation  in  the  nose  region  in  the  second  flight  behaved  in 
the  same  as  manner  as  in  the  First  flight,  though  at  a  different  level. 

Significant  radiation  was  detected  both  from  the  centerline  and  the  far-field 
regions  of  the  plume. 


B.  Progress  of  the  Work 

Before  the  first  flight  it  has  been  predicted,  based  on  existing  experimental 
and  theoretical  d?‘.  that  the  flowfield  will  be  in  thermochemical 
nonequilibrium.  ‘Iter  the  first  flight  experiment,  calculations  were  carried 
out  for  the  radiative  signal  received  in  the  blunt  nose  region  of  the  vehicle 
using  the  then  standard  nonequilibrium  codes.  The  standard  codes  consisted  of 
the  t'.o-temperature  nonequilibrium  blunt  body  flow  code  and  the 
nonequilibrium  radiation  code  (NEQA1R).  In  the  low  altitude  range,  namely, 

1  etween  40  and  50  km,  the  calculated  radiation  intensity  agreed  approximately 
with  the  flight  data.  However,  at  higher  altitudes,  the  calculation 
underpredicted  the  measurement  by  at  least  seven  orders  of  magnitude.  This 
trend  was  later  found  to  be  the  same  even  for  the  second  flight. 

This  discrepancy  between  the  measured  and  the  calculated  radiation  intensity 
of  NO  set  off  concerted  efforts  within  the  Branch,  IDA,  and  elsewhere.  The 
immediate  response  was  to  attempt  explaining  the  experimental  data  in  terms 
of  equilibrium,  that  is,  that  the  flowfield  and  radiation  mechanism  were  in 
thermal  equilibrium.  The  equilibrium  calculation  yielded  radiation  intensity 
values  that  are  uniformly  an  order  of  magnitude  higher  than  the  flight  data  at 
all  altitudes.  The  fact  that  this  discrepancy  existed  at  altitudes  of  40  to  50  km, 
where  NEQAIR  gave  a  satisfactory  result,  was  particularly  disturbing.  This 
reinforced  the  belief  that  the  flowfield  is  in  nonequilibrium. 

It  was  understood  at  once  that  the  original  version  of  the  code  NEQAIR  is 
inappropriate  for  this  problem  because  the  intended  flight  speeds  are 
different:  the  code  assumes  that  electrons,  prevelant  only  at  high  flight  speeds 
(above  5  km/s),  are  responsible  for  excitation  of  internal  states.  Levin  made  a 
correction  to  the  NEQAIR  code  to  effect  that  it  is  the  heavy  particles,  and  not 

electrons,  that  cause  excitation  of  electronic  states.  This  correction  raised  the 

calculated  radiation  intensity  at  high  altitudes  by  an  order  of  magnitude, 
leaving  the  discrepancy  to  be  two  orders. 

Thereupon,  several  hypotheses  were  developed  to  explain  the  still  remaining 
discrepancy  between  the  theory  and  the  flight  data.  These  included:  (1) 
neglecting  of  Burnett  terms  in  the  constituency  equations  in  fluid  mechanical 
equations,  (2)  slip  at  the  wall  boundary,  (3)  inaccuracy  in  modelling  of  the 

rotational  relaxation  behavior,  (4)  inaccuracy  in  the  value  of  the  cross  section 

for  excitation  of  electronic  states  of  NO,  (5)  inaccuracy  in  the  spectroscopic 
and  radiative  parameters  for  NO  band  systems,  (6)  inaccuracy  in  modelling  the 
NO  formation  processes,  and  (7)  the  nonBoltzmann  distribution  of  translational 
mode  in  the  shock  wave. 

Lumpkin  et  al  first  explored  the  effect  of  neglecting  of  Burnett  terms.  They 
found  that  inclusion  of  Burnett  terms  improved  the  theory  slightly,  that  is,  by 
an  order  of  a  factor  of  1.5. 

Lumpkin  also  pioneered  modelling  of  the  rotational  relaxation  correctly.  This 
was  done  by  extending  the  existing  rotational  model  by  incorporating  the 
recent  theoretical  and  experimental  data  obtained  from  the  linewidth 
measurements.  This  correction  produced  a  slight  improvement.  Later,  Moreau 
et  al  made  a  drastic  modification  to  rotational  model  by  incorporating  the  so- 
called  diffusion  model.  This  correction  greatly  increased  the  calculated 
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radiation  intensity. 

The  cross  section  for  electronic  excitation  of  NO  molecule  was  studied 
experimentally  at  Stanford  Research  International  (SRI).  The  work  showed 
that  the  cross  section  is  very  large,  that  is,  about  50  square  angstroms. 

Boyd  et  al  explored  the  question  of  whether  the  non-Boltzmann  energy 
distribution  is  responsible  for  the  discrepancy  between  the  flight  data  and  the 
continuum  calculations  using  the  Direct  Simulation  Monte  Carlo  (DSMC) 
technique.  It  was  shown  first  that  DSMC  technique  is  at  least  equivalent  to  the 
continuum  approach  in  the  regime  where  they  are  expected  to  be,  and  that 
indeed  the  non-Boltzmann  distribution  is  partly  responsible  for  the 
discrepancy.  But  the  extent  of  this  effect  was  found  to  be  uncertain. 

Laux  et  al  and  Moreau  and  Laux  investigated  the  accuracy  of  the  spectroscopic 
and  radiation  intensity  parameters  for  NO  by  comparing  the  calcualted  results 
with  the  experiment  conducted  at  Stanford  University  with  a  plasma  torch.  By 
using  the  spectroscopic  and  radiative  parameters  generated  recently  at 
Computational  Chemistry  Branch,  rather  than  the  set  used  in  the  original 
version  of  NEQAIR,  they  were  able  to  demonstrate  that  the  experimental  data 
can  be  closer  reproduced.  This  led  to  creation  of  the  NEQAIR-2  code.  Use  of  this 
code  narrowed  the  discrepancy  still  further. 

Moreau  developed  a  five-temperature  model  to  describe  the  process  of 

formation  and  removal  of  NO  molecules  in  an  experiment  conducted  at 
CALSPAN  (W.  Wurster's  experiment).  The  model  so  developed  can  numerically 
recreate  CALSPAN  experiment  by  assigning  a  vibrational  temperature  of  NO 
molecule  which  is  different  from  that  of  N2  molecule.  The  new  model  lowers 
the  calculated  radiation  intensity  at  high  altitudes. 

The  last  two  innovations  by  Moreau  have  not  yet  been  fully  incorporated  in 
the  calculation  of  radiation  intensity.  However,  if  one  assumes  that  the  cross- 

section  of  excitation  of  NO  electronic  states  is  infinitely  large,  then,  by 
incorporating  all  the  improvements  made,  the  calculated  intensity  agrees 

closely  with  the  flight  data  at  high  altitudes,  for  the  first  flight.  For  the  second 
flight,  the  same  approaches  overestimates  the  radiation  intensity. 

In  order  to  study  the  dynamics  of  the  plume,  a  theoretical  calculation  was  first 
made  of  the  fluid  flow  in  the  plume.  The  flight  data  on  the  plume  region  were 
calculated  and  compared  by  Candler  el  al.  A  problem  exists  in  that  the 
thermochemical  state  of  the  rocket  effluent  at  the  nozzle  exit  is  uncertain.  By 
adjusting  the  state  parameters  at  the  exit,  it  was  shown  that  the  observed 
radiation  intensity  along  the  centerline  can  be  numerically  recreated.  The 

radiation  is  attributed  to  that  from  the  solid  particulates  in  the  effluent. 

However,  for  the  far-field  radiometer,  that  is,  that  viewing  at  25  degrees,  the 
calculation  severely  underpredicted  the  measured  value.  The  experimental 
data  can  be  recreated  numerically  by  first  assuming  that  the  radiation 
emanates  mostly  from  the  CO  Cameron  band  system,  and  that  the  relative 

populations  of  electronic  states  of  CO  molecules  are  frozen  at  the  nozzle  exit 

conditions.  Thus,  the  far-field  radiation  signal  seems  to  emanate  from  the  gas- 
phase  radiation. 
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C.  Summary  of  Most  Important  Results 

1.  It  was  shown  that,  by  including  the  Burnett  terms  in  the  constituent 
equation,  accounting  for  the  non-Boltzmann  energy  distribution,  correctly 
modelling  the  NO  formation  and  rotational  excitation  processes,  and  using  the 
correct  transition  probabilities  for  NO  radiation,  one  can  numerically 
reproduce  the  observed  radiation  of  NO  in  the  bow  shock  region  for  the  flight 
velocity  of  3.5  km/s,  provided  the  rate  of  excitation  of  the  electronic  states 
through  the  vibration-electronic  energy  transfer  is  very  fast.  (For  the  5  km/s 
flight,  the  agreement  between  calculation  and  measurement  is  not  as  good.) 

2.  It  was  found  that  the  DSMC  technique  is  at  least  equivalent  to  the  continuum 
technique  in  determination  of  nonequilibrium  radiation. 

3.  In  the  expanding  plume  flow  behind  the  second  flight  vehicle,  the 
electronic  excitation  temperature  seems  to  be  frozen  at  the  value  at  the  exit  of 
the  nozzle. 
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